Multipartite continuous- variable entanglement from concurrent nonlinearities 
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We show theoretically that concurrent interactions in a second-order nonlinear medium placed 
inside an optical resonator can generate multipartite entanglement between the resonator modes. We 
show that there is a mathematical connexion between this system and van Loock and Braunstein's 
proposal for entangling N continuous quantum optical variables by interfering the outputs of TV 
degenerate optical parametric amplifiers (OPA) at a TV-port beam splitter. Our configuration, 
however, requires only one nondegenerate OPA and no interferometer. In a preliminary experimental 
study, we observe the concurrence of the appropriate interactions in periodically poled RbTiOAsCU. 

PACS numbers: 03.67.Mn, 03.65.Ud, 03.67.-a, 42.50.Dv, 42.65.Yj 



A crucial direction of effort in quantum information 
is the entanglement of (many) more than two systems. 
The current record number of entangled q-bits is four 
0. Continuous variables (CV) are a fascinating alterna- 
tive to discrete ones and lend themselves well to quan- 
tum optical implementation |2j. Multipartite entangle- 
ment of continuous variables was proposed by van Loock 
and Braunstein 0, Q and experimentally demonstated 
in two different regimes d E|. A CV multipartite en- 
tangled state is an inseparable multimode squeezed state 
that tends towards a Greenberger-Horne-Zeilinger (GHZ) 
state in the limit of infinite squeezing. In van Loock 
and Braunstein's method, such states are generated from 
N squeezed modes of the field emitted by optical para- 
metric oscillators (OPO's) below threshold [i.e. optical 
parametric amplifiers (OPA's)] and combined by appro- 
priately balanced beam splitters. The quantum interfer- 
ence at the beam splitters requires interferometric stabi- 
lization of the optical paths and indistinguishability, i.e. 
frequency and polarization degeneracy, of all N modes. 
In this paper we show that multipartite entanglement is 
obtainable by use of a single OPA and no beam split- 
ters. The OPA nonlinear medium must simultaneously 
phase-match several second-order nonlinearities that cre- 
ate two- mode squeezing between N cavity field modes. 
The advantage of this scheme for experimental purposes 
is that it can be made very compact in a single period- 
ically poled ferroelectric nonlinear crystal, with no in- 
terferometer to lock. Moreover, there is no degeneracy 
constraint on the frequencies. In the next section of this 
paper, we expose the theoretical arguments that prove 
our assertion. We then detail how a concurrence of three 
nonlinearities, suitable for entangling four modes, can be 
created in periodically poled ferroelectrics and we present 
preliminary observations of such a concurrence in period- 
ically poled RbTiOAs0 4 (PPRTA). 

It is well known that a two-mode squeezer such as a 
nondegenerate OPA can generate an EPR state out of 
vacuum or coherent input 0,0- The Hamiltonian in the 
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interaction picture is 

H\ = ih(3x{a\a\ — 0102), (1) 

where x is t ne nonlinear coupling coefficient and (3 the 
real, assumed undepleted, coherent pump field ampli- 
tude. Solving the Heisenberg equations for the fields 
gives the squeezed joint quadratures 

Pi(t)+P 2 (t) = (Pi+P 2 )e~^ t (2) 
Xi(t)-X 2 (t) = (Xi-X 2 )e-^ t (3) 

where X = (a + a^)/V2, P = i(at - a)/V2, and 
a(t = 0) = a. These EPR operators commute and admit 
maximally entangled common eigenstates such as 

/ \x)i\x) 2 dx= / \p)i\ -p) 2 dp= 5^|n)i|n) a (4) 

^ "* 71 = 

in the limit of infinite squeezing 0xt — > 00 0, . We 
now ask whether concurrent interactions involving three 
modes would yield a tripartite CV entangled state, e.g. 
J \xxx)dx. It is simple to check that the Hamiltonian 

H 2 = ih(3xa{a\a\ - aia 2 ) + ihf3xb{a\a\ - a 2 a 3 ), (5) 

does not create such CV tripartite entanglement: the so- 
lutions of the Heisenberg system yield only two squeezed 
joint modes out of three, the third one being a constant 
of motion initially subject to vacuum fluctuations. More 
interesting is the symmetrized three-mode Hamiltonian 
(we now take the interaction strengths equal, for the sake 
of simplicity and symmetry) 

H 3 = ihPx{a\a\ + a\a\ + 4<4) + H.c, (6) 

whose system of Heisenberg equations is A = M.A', 
where A T = (01,02,03). M. has only zeroes on 
the diagonal, j3\ everywhere else, and eigenvalues 
(2/3%, — 0x> ~Px)- The eigenmodes are joint operators 

Pi(t) + P 2 (t) + P 3 (t) = (Pi + P 2 + P 3 )e- 2 ^ (7) 
Xi{t)-X 2 {t) - (Xi - X 2 )e-^ 1 (8) 
Xi(t)-X 3 (t) = (Xi-X 3 )e-^\ (9) 
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whose common eigenstate is a multipartite entangled 
state that tends towards the GHZ state J \xxx)dx when 
fixt — ► oo. It is straightforward to generalize to 



N N 



H N = ih/3 x a H + i/ ' c -' 



(10) 



=1 j>i 



which gives the set of multipartite entangled modes 



N 



N 



i=l 

X|(t)-X,(t) 



(11) 



(Xi-X&Vi^j. (12) 



(Note that the phase sum squeezing is higher than for 
the other GHZ modes.) It would therefore seem that we 
have found a procedure for entangling an arbitrary num- 
ber of modes: one just needs to have every one of them 
interacting equally strongly with all the others. Before 
we turn to the experimental feasibility of this scheme, we 
examine its relation to that of van Loock and Braunstein. 

A connexion is already known for bipartite entan- 
glement. Indeed, Hi [Eq. QJ] can also be viewed as 
the result of the transformation of the Hamiltonian of 
two single-mode squeezers, H[ = 'i^Px[~{ a i < — a i) + 
(a^ — a!)], by a lossless balanced beam splitter: H'{ — 
UssHiUgg = Hi Both possibilities have been used 
experimentally: Hi in the first CV EPR experiment 
and subsequent ones 1 1 2L 1 1 31 and H" in quantum tele- 
port at ion experiments [lj, [l5j . 

We now show the same kind of simple relationship 
exists for CV multipartite entangled beams. The orig- 
inal proposal of CV multipartite entanglement uses 
three single-mode squeezed beams, produced by H' 3 = 

*f 0X(~ a i + a V + a 3 t ) + H.c. and mixed by a "tritter", 
i.e. the combination of a 2:1 and a 1:1 beamsplitter 
The transformation of H$ by the tritter is 



HZ = - 

3 3 



(2t 



b$) + H.c. - 2H 3 



(13) 



i.e. a combination of our H3 [Eq. ©] and symmetrized 
single-mode squeezers. The two-mode dependence of HZ, 
and H3 is therefore identical. 

The corresponding TV-mode entangling Hamiltonian is 



N-2 
2N 



N 



H.c. - -H N , 



(14) 



where Hm is our entangling Hamiltonian [Eq. (| X(jp] . The 
system of Heisenberg equations for H'^ is B = M.B', 
where B T = (pi, ... , &zv). All of A^'s diagonal elements 
equal (N — 2)P X /N and all off-diagonal elements equal 
—2(3 X /N. The eigenvalues are (—flx, 0X> ■ ■ ■ > fix) an d 
eigenvectors of the form of Eqs. 1)11112(1 . with X and 
P swapped and equal squeezing rates exp(— /?xi). Our 
Hamiltonian Hn (or ^3) leads to the same matrix, but 
with zero diagonal, and asymmetric squeezing rates [Eqs. 



(|llll2fl ] . The advantage is that no control is required over 
the phases of the input squeezed states, i.e. over the signs 
of the different terms in H[ 3 N . This greatly simplifies 
the experimental setup, as one passes from the interfer- 
ence of N OPA's to the output of a single OPA. 

The experimental principle is to use the simultaneous 
nonlinear interaction of different eigenmodes of an opti- 
cal resonator. The presence of the resonator here simpli- 
fies the situation by selecting a discrete comb of resonant 
modes out of the quantum field continuum. Moreover, 
because the assembly of the optical resonator with the 
nonlinear medium constitutes a parametric oscillator, it 
provides us with the option of operating either in the 
spontaneous emission regime (below oscillation thresh- 
old) or in the stimulated emission regime (above oscilla- 
tion threshold). We plan to investigate the latter case in 
subsequent work but we restrict the scope of this paper 
to vacuum-seeded optical parametric amplification. 

In general, a nonlinear second-order medium pumped 
at frequency lo p will optimally couple pairs of modes 0^2 
such that u>i -\- u>2 = w p (from the phase-matching con- 
dition). The phase-matching bandwidth can be broad 
enough (e.g. 100 GHz) for several pairs of modes, sepa- 
rated by a free spectral range (FSR), say, A < 1 GHz, 
to have approximately the same coupling strength. It is 
also well known that the resonance condition depends on 
the frequency, via dispersion, but these effects are small 
enough to be negligible over a few free spectral ranges. 
Still, a singly-pumped OPA cannot realize multipartite 
entanglement. The two possible cases are: 

(i) uj p coincides with twice a cavity resonance frequency 
2bj(j. Then the Hamiltonian is 



,t„t 



+ aja_i + a!>a_-) + 



t„t 



2"-2 



H.c. 



(15) 



where ao has frequency u>q and a±k, to±k = <^o + ±fcA. 

(ii) w' coincides with the sum of two consecutive cavity 
resonance frequencies, e.g. lu' q = u>q + u>i. The Hamilto- 
nian is 



ihPx (^ a o a i 



+ ali at + a 



-2 U 3 



H.c. 



(16) 



In neither case is a multipartite entangling Hamilto- 
nian realized because we never have a given set of more 
than two modes all connected together by the interac- 
tion. This problem can be easily solved by having sev- 
eral pump beams at different frequencies, more precisely 
half a FSR apart. This amounts to adding both Hamil- 
tonians l|15fl and i|16fl . However, an additional difficulty 
stems from the fact that the degenerate interaction in 
(|15fl has the same sign as the other nondegenerate inter- 
actions, contrary to what is required by Eq. I|14|). In fact, 
it is impossible to have a different interaction phase for 
two different downconversion terms that share the same 
pump such as in 



ihPx 



H.c. 



(17) 
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to be compared with Eq. I|15fl • This makes the exact re- 
alization of van Loock and Braunstcin's Hamiltonian H'^ 
impossible with concurrent interactions, as it requires op- 
posite signs for the degenerate and nondegenerate inter- 
actions [Eq. (|14() ]. Thus, one must go back to the other 
side of the multiport splitter, where the individually tun- 
able optical paths provide the necessary degrees of free- 
dom. However, concurrent interactions can realize our 
alternate, nondegenerate Hamiltonian [Eq. (|10fl ]. 

In practice, the implementation of requires obtain- 
ing all possible nondegenerate coupling terms without 
any degenerate one, for a given set of modes. We now 
show how this can be achieved for sets of three or four 
modes. Several equivalent possibilities exist. Without 
tediously enumerating them all, we focus, without loss 
of generality, on the simplest ones that are experimen- 
tally realizable. These are the only ways we have found 
to fulfill the two requirements stated above (which does 
not, obviously, constitute a proof of unicity) . We assume 
propagation along the principal axis £ of a nonlinear crys- 
tal and use the polarization (y, z) and frequency degrees 
of freedom to label the modes. Graphical representations 
of the interactions are given on Fig. ^ Three-mode en- 




FIG. 1: concurrent entangling interactions (dashed lines) for 
(a): three and (b): four OPA modes (top arrows). Bottom 
arrows represent the pump fields. 

tanglement is realized by simultaneously phase-matching 
zzz and yzy parametric downconversion [Fig. ^a)]: 

H 3 = ih [/3 y (2uJo)xyzyal(uJo)ai(ujo) 

+ Py(w + Wl)Xy*3/<4(Wo)4( W l) 

+ f3 z (uj + wi)x*«ol(wo)ol(wi)] + H.c. (18) 

No other interaction occurs inside the set of modes 
{a v (uJo), a z (wo), a z (a>i)} [in particular, yyy is not phase- 
matched and there is no /3 z (2wo)]- The pump ampli- 
tudes can be adjusted to make up for residual differ- 
ences in the nonlinear coefficients Xijk so that the inter- 
action strengths are all equal. Under these conditions, 
this nondegenerate OPA should amplify vacuum inputs 
into tripartite entangled modes of the form Q. The 
four-field entanglement Hamiltonian necessitates ad- 
ditional phase- matching of the yyy interaction, which can 



be done by using a zzz crystal rotated by 90°. Neglecting 
the FSR difference between the two polarizations (which, 
for practical purposes, can be made arbitrarily small us- 
ing two crystals) we can design an entangling interaction 
between four equally spaced cavity modes [Fig.QIb)]: 

Hi =iK{(3y{u Q + ui)xyzya\{uo)a\{iJi) 
+ (3y(2u)i)xyy y ay(tv )al(uJ2) 

+ f3 z (2uj2)Xzzzai(ui)al(uj 3 ) 

+ (3y{u 2 + uj 3 )xyz y al(uj2)al(uj 3 )} +H.c. (19) 

Again, these are all the possible pair couplings inside the 
set of modes {a y (uJo), a z {^>\), a y (uj2), a z (^ 3 )} and these 
should thus all be entangled by H4. Note that the ab- 
sence of yzz and zyy prevents degenerate terms from 
appearing. The equidistant pump frequencies can eas- 
ily be obtained by acousto-optic or electro-optic modula- 
tion. This scheme does not scale easily to N > 4, unfor- 
tunately, because this requires phase-matching zyy /yzz 
and it then becomes impossible to avoid degenerate in- 
teractions with the wrong sign, which, we have found, al- 
ways dramatically reduce the number of entangled modes 
(detailed calculations will be published elsewhere). This 
is a consequence of the freezing of the optical phases 
in our scheme; however, this very limitation is precisely 
what makes possible the experimental simplicity of our 
approach, which could lead to extremely compact and 
low-loss CV multipartite entanglers. 

As we have seen, realizing H 3 ^ in a nonlinear opti- 
cal material necessitates the concurrence of two different 
interactions. Using two different nonlinear crystals is a 
possibility but it is costly in terms of optical losses, espe- 
cially inside an optical resonator. It is therefore desirable, 
in line with the rationale of this paper, to obtain all in- 
teractions in the same crystal. While such coincidences 
are extremely rare in birefringent phase-matching, they 
are very easy to engineer in quasi-phase-matched materi- 
als. Quasi-phase-matching (QPM) is as old as nonlinear 
optics It relies on spatial modulation of the nonlin- 
ear coefficient, here by periodically poling a ferroelectric 
crystal, to make up for the phase mismatch of a particular 
nonlinear interaction 01 • QPM allows practically any 
interaction to be phase-matched in the same material by 
simply changing the poling period. Poling the same crys- 
tal with two different periods therefore yields the desired 
coincidence. Even simpler designs are possible: one can 
find two interactions sharing the same period |l8j| . More- 
over, since the spatial modulation of the nonlinear coeffi- 
cient is a square wave, one can also find coincidences be- 
tween different poling harmonics. We have observed one 
instance of these in PPRTA. RTA is an isomer of KTP 
with lower residual absorption losses and equivalent non- 
linear coefficients. Like KTP, RTA is non-hygroscopic, 
has a very high optical damage threshold, and presents 
neither photorefractive damage nor blue- induced infrared 
absorption. We consider the particular set of Nd-doped 
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FIG. 2: Simultaneous quasi-phase-matching of yzy and zzz 
SHG in PPRTA. The SHG powers are measured versus crys- 
tal temperature for the zzz and yzy interactions. Concur- 
rences involving secondary maxima are at 18.75°C and 57.5°C 
(dashed lines). Inset: full-scale data plot. 

laser wavelengths of 532 nm, for the pump fields, and 
1064 nm, for the parametrically amplified fields. The pol- 
ing periods required to quasi-phase-match SHG of 1064 
nm at room temperature are 43 /im for yzy and 8.37 /mi 
for zzz. The interaction is tunable by varying the re- 
fractive indices via the wavelength, incidence angle, or 
temperature. The goal is to obtain simultaneous QPM 
of both yzy and zzz interactions at the same tempera- 
ture. The 8.37 /im period is the fifth harmonic of 41.85 
/im, close to 43 /im. We must therefore find an inter- 



mediate poling period for which both interactions occur 
at the same temperature. The temperature dependance 
of RTA indices being not perfectly known yet, we have 
tried a 41.95 /im crystal. The results of SHG measure- 
ments are plotted in Fig. [21 The input beam at 1064 nm 
and the output beam at 532 nm are respectively polar- 
ized and analyzed with polarizers. Even though the two 
main QPM peaks are still separated by about 40°C, we 
already have coincidences of each main peak of one in- 
teraction with a secondary one of the other interaction. 
The ratio of the two maxima gives {dy^ y /d e Jl z ) 2 = 7, for a 
theoretical value of (Sdzi/dzz) 2 ~ 1.8. We attribute the 
discrepancy to high frequency irregularities in the crys- 
tal's poling, which would affect the fifth harmonic of the 
period more than its fundamental. We are currently re- 
fining measurements of temperature tuning for RTA in 
order to improve the period design. 

In conclusion, we have demonstrated that concurrent 
nonlinearities in a single OPA can entangle four modes, 
possibly more. The current technology makes such an 
OPA quite feasible and we are now preparing an exper- 
imental realization. The compactness and simplicity of 
such a source of entangled beams should make it very at- 
tractive for CV quantum information implementations, 
such as teleportation networks, controlled dense coding, 
and quantum error correction via telecloning. 

We thank Richard Stolzenberger and Tito Ayguavives, 
both formerly of Coherent Crystal Associates, Inc., for 
their help with PPRTA. This work is supported by NSF 
grants PHY-0245032 and EIA-0323623. 

Note added in proof: we have become aware of related 
work for tripartite photon- number entanglement |l9j| . 
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